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Abstract. It is well known that long-range interactions pose serious problems for the formulation of statis- 
tical mechanics. We show in this paper that ensemble equivalence is violated in a simple mean-field model 
of A'' fully coupled classical rotators with repulsive interaction (antiferromagnetic XY model). While in the 
canonical ensemble the rotators are randomly dispersed over all angles, in the microcanonical ensemble a 
bi-cluster of rotators separated by angle n, forms in the low energy limit. We attribute this behavior to 
the extreme degeneracy of the ground state: only one harmonic mode is present, together with TV — 1 zero 
modes. We obtain empirically an analytical formula for the probability density function for the angle made 
by the rotator, which compares extremely well with numerical data and should become exact in the zero 
energy limit. At low energy, in the presence of the bi-cluster, an extensive amount of energy is located in 
the single harmonic mode, with the result that the energy temperature relation is modified. Although still 
linear, T = aU, it has the slope a ~ 1.3, instead of the canonical value a = 2. 

PACS. 05.20.-y, Classical statistical mechanics - 05.20.Gg Classical ensemble theory - 05. 45. -a Nonlinear 
dynamics and nonlinear dynamical systems 



1 Introduction 

The relation between microscopic dynamics and macro- 
scopic thermodynamic behaviour can nowadays be explored 
in full detail in computer simulations Q. This allows one 
to test important hypotheses of statistical mechanics. For- 
most among these is equivalence of ensembles. 

It is widely accepted that the constant energy micro- 
canonical ensemble gives the same results for average val- 
ues as the constant temperature canonical ensemble. This 
is true under certain conditions |^ , among them the most 
important is that interactions must be short-ranged. If 
interactions are long-range and attractive, as for gravitat- 
ing systems, all of thermodynamics breaks down due to 
the non extensivity of thermodynamic potentials For 
systems of Coulomb charges of opposite signs, screening 
effects may help in the construction of thermodynamics, 
but the problem is not trivially solved, even for classical 
systems and each case must be examined separately. 

In this context, mean- field models occupy a special sta- 
tus. Here, thermodynamic potentials can be made exten- 
sive if the thermodynamic limit is performed by rescal- 
ing the coupling with system size and letting the range 
of the interaction go to infinity j|] . Mean- field models are 
quite well studied in the canonical ensemble and serve as a 
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zero-order characterization of phase transitions. Although 
their solution is often trivial for systems without disor- 
der, they may hide important subtleties for disordered 
and frustrated systems . Far fewer studies of mean-field 
models exist in the microcanonical ensemble. This may be 
the reason why, until now, ensemble equivalence has not 
been questioned for these models. In fact, one might have 
already expected some surprises on the basis of the exact 
solution of a model by Hertel and Thirring ||^, where in 
the mean-field limit in the presence of extensive thermo- 
dynamic potentials ensemble inequivalence was explicitly 
shown. 

In this paper we present a detailed study of the low 
temperature/energy phase of a model of classical rotators, 
whose potential energy is that of the mean-field antifer- 
romagnetic XY model. The ground state of this model 
is highly degenerate, and while in the canonical ensem- 
ble equilibrium states are disordered at all temperatures 
(i.e. rotators do not display any directional organization), 
microcanonical ensemble simulations show the presence 
of a "bimodal" state where rotators are mainly grouped 
in two "clusters", pointing in directions separated by an 
angle tt. This dynamical effect also has thermodynamical 
consequences: the order parameter measuring the degree 
of clustering is non-zero in the limit of zero energy and 
the energy temperature relation is not that predicted by 
the canonical ensemble. We find that, although the energy 
temperature relation T = all is linear, as in the canonical 
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case, the coefficient is a w 1.3 and not the canonical value 
a = 2. From the point of view of dynamics, the rotators 
can be separated into two groups: a slow group which os- 
cillates around the bi-cluster and a group of almost freely 
rotating rotators, which we call a "gas" , following an anal- 
ogy with particle motion. 

In Section ^ we introduce the model. In Section ^ we 
discuss the unusual properties of its ground state and in 
the following Section ^ we present the main controversial 
points related to ensemble inequivalence. In Section | we 
present an analytical expression for the probability den- 
sity function (PDF) for the orientation of a rotator, while 
the system is in the bi-cluster state, and we explain the 
consequences for the statistical properties in the system in 
the microcanonical ensemble. In Section ^ we discuss the 
main features of the dynamics of the moments of the PDF. 
The paper ends with some conclusions and perspectives. 



2 The model 

We consider classical rotators denoted by the angle 9i, 
i — I, . . . , N, which all interact with each other, with an 
antiferromagnetic coupling and an external field h 



(1) 



^,J 



We will mostly restrict ourselves to h — 0. 

After defining the complex order parameters (with i 
the imaginary unit) 

Mfe = ^ ^ expitkOn) = \Mk\ exp(iVfc) , (2) 



the potential can be rewritten as 

V (^-^~h\M,\sin^P,^ (3) 

The microcanonical ensemble is obtained by adding a ki- 
netic energy term to the above potential 



H = K + V 



with 



K 



E Pn 
2 

n=l 



(4) 



(5) 



In this formulation the model can also be thought of as de- 
scribing a system of particles with unitary mass, interact- 
ing through the mean-field coupling V (the names rotator 
and particle will be equivalently used throughout the pa- 
per). This model has been first introduced in Ref. |^ for 
2D nearest-neighbour couplings and then studied in the 
antiferromagnetic mean-field context in Ref. 

The total energy E = UN {U being the energy den- 
sity) is fixed by the initial conditions and is conserved 
in time, while temperature is defined through the time 



averaged kinetic energy (see Ref. ||T^ for alternative defi- 
nitions of temperature), T = 2 < K > /N, where < • >= 



lim - / . If /i = 0, as is usually the case, then the total 
momentum P = p„ is also conserved. We set P = in 

n 

order to avoid ballistic center of mass motion when h ~ Q. 
The equations of motion 



\Mi \ [sin(0„ — ■i/'i) -f hcos9„ 



(6) 



have been integrated using an improved fourth-order sym- 
plectic scheme The algorithm is 0{N), provided one 
first computes Mi in the central loop. During the time 
evolution, we sample the instantaneous values of the order 
parameters @) up to fc = 20 and we compute their running 
time averages. System size was varied from iV = 100 to 
N = IQ^. We have also performed canonical Monte-Carlo 
simulations, using the Metropolis algorithm, for compari- 
son. 

The initial conditions were of two classes: i) homoge- 
neous state On = {2TTn)/N, which can be shown to be 
marginally stable (see below), to which we add either a 
small spatially random perturbation 9n + Vn and/or 
a small momentum p„ = r„ with zero average; ii) homo- 
geneous state On = {2Tm)/N with p„ = A sin 9m which 
leads to a faster induction of the bimodal state we want 
to study. These initial conditions lead to the same bimodal 
state, discussed below, at low energy. 

The canonical solution of this model for /i = is 
sketched in Ref. It uses the Hubbard- Stratonovich 
trick to decouple the rotators and the thermodynamic 
limit is performed by a saddle-point technique. The re- 
sult is that all moments (^) vanish in the N ^ oo limit, 
including the magnetization A/i, which implies, on the one 
hand, that 

T = 2U , (7) 



on the other, that since in this limit 



Ah 



V{0) exp{ik9)d9 



(8) 



the PDF V{9) is flat. Hence, the bimodal state is absent 
in the canonical ensemble, which we have confirmed by 
Monte-Carlo simulations. 

Most of the discussions below will concentrate on the 
reasons for the different findings in the microcanonical and 
canonical statistical ensembles. 



3 Ground State and Statistical Properties in 
the Canonical Ensemble 

The long ranged interactions mean that it is impossible 
to satisfy all the antiferromagnetic bonds at once and the 
model is highly frustrated. In zero field, the frustration 
is minimized for configurations with Mi — 0, giving a 
ground state energy of [/ = E jN = 0. 
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The ground state is infinitely degenerate, as there is 
an infinity of ways of minimizing the frustration. For ex- 
ample, grouping the rotators into pairs with angles 9i and 
Oi + TT ensures that Mi = for all configurations of the 
pairs. The pairs do not have to be arranged in an ordered 
way however, as any disordered arrangement will equally 
give Ml = 0. Neither is the ground state manifold re- 
stricted to pairs: one can equally construct ground states 
from groups of three, four, five ... rotators separated by 
angles of 27r/3, 27r/4, 27r/5.... By moving rotators in clus- 
ters whose total angle is zero the system can evolve from 
one ground state to another remaining on the constant 
energy hypersurface. 

The high dimensional ground state manifold follows 
from the fact that the ground state condition requires the 
two constraints AIi^ ~ Miy — only. One can there- 
fore expect a ground state to possess N — 2 unconstrained 
degrees of freedom, which is easily confirmed by calcu- 
lating the Hessian Jij — —d^V/dOtdOj. For example, for 
the perfectly homogeneous ground state at ft. = 0, with 
6*, = {2m) /N 



J, 1 = cos 

2iV 



Mi - j) 
N 



(9) 



The matrix indeed has two non-zero equal eigenvalues, 
— 1/4 and N — 2 zero eigenvalues. 

Collective organization of the particles into reduced 
symmetry states can reduce the number of constraints 
even further. For example, in our system, if the rotational 
symmetry is broken and the spins lie along a single axis 
the ground state condition is reduced to a single constraint 
= 0, leading, at the harmonic level to TV — 1 free de- 



grees of freedom W2\ 
lating the Hessian 



This result is confirmed by calcu- 
J) , on a perfect bi-cluster groundstate 
with N/2 rotators at 6* = and N/2 at 6* = tt. One finds 
only one non-zero eigenvalue, —1/2, corresponding to the 
counter vibrating motion on the circle of the two groups 
of particles; all the other eigenvalues are zero. 

The model we study is an extreme case of a collective 
paramagnet [|l^, or classical spin liquid [ p^|p^ ; a system 
that remains disordered with no evidence of spin freez- 
ing down to the limit of zero temperature. The special 
points on the ground state manifold with a reduced num- 
ber of constraints can dominate the partition function, in 
the canonical ensemble, leading to an "Order by Disorder" 
transition Jisf to a reduced symmetry state. However, the 
mode counting arguments of Moessner and Chalker |^ 
predict this to happen only if the number of liberated 
modes at the special points exceeds the number of zero- 
modes in a state with full symmetry. This is certainly 
not the case here and as we have confirmed by Monte 
Carlo simulation. Order by Disorder does not occur. We 
have simulated between TV = 10 and N = 10^ rotators 
down to temperatures T < 10~^. At all temperatures, the 
system remains perfectly disordered, with no evidence of 
bi-cluster formation. The number of zero modes can be di- 
rectly verified by measuring the specific heat at constant 
field, C/i, at low temperature, as each quadratic mode 
makes a contribution 1/2 (in units of ks), while each zero 



mode makes a contribution zero. We find Ch = 1.0, for 
the N rotator system, as expected for two regular modes. 
There is therefore no evidence of the system preferring 
states with a single quadratic mode. 



4 Statistical Properties in the Microcanonical 
Ensemble 

The system reserves a surprise, when studied in the micro- 
canonical ensemble, as we do not observe ensemble equiv- 
alence. For the classes of initial conditions cited above, 
the low energy state of the system is not one with a ho- 
mogeneous distribution of angles, rather we observe the 
formation of a bimodal structure, with enhanced prob- 
ability for two angles separated by a distance tt As 
the energy goes to zero, the asymptotically reached state 
has broken symmetry, but is not perfectly bi-modal, as 
we discuss below. Nevertheless, when the Hessian for such 
a state is calculated numerically, we find that it also has 
only one non-zero eigenvalue. It seems therefore that the 
formation of the bicluster in the microcanonical ensem- 
ble, realises the condition of minimizing the number of 
non-zero modes. 

After a transient time r, the formation of a stable 
bi-cluster is revealed by a non-zero value of the second 
moment IM2I of 'P{9). The early time evolution of IM2I 
shows the typical behavior observed for the growth of the 
mean-field in self-consistent dynamical models |]l6| with 
an initial exponential growth, followed by a saturation 
reached after a darnped oscillatory motion has died away, 
see Figs. |l|(a) and |^(b). In Figs. |l|(c) and |l|(d) we show 
the space-time plot of V{9,t), which reveals the origin 
of the oscillations in the density waves originated peri- 
odically from the bi-cluster, preceding stabilization. The 
darkest colours correspond to the highest densities. A very 
sharply peaked but unstable bi-cluster forms over a short 
time period. The structure disperses with well defined in- 
stability edges that propagate from the cluster center and 
re-appears in a quasi-periodic manner with slowly length- 
ening period. The dispersing cluster appears to interact 
with propagating fronts from previous incarnations. The 
result is that the dispersion is successively slower and the 
fronts less well defined for following quasi-periods until fi- 
nally the bi-cluster stabilizes. In addition Fig. [|(d) shows 
that once this coherent structures has emerged and is sta- 
bilized, it propagates around the circle with apparently 
ballistic dynamics. 

In Fig. ||(a), we show IM2I, averaged over times later 
than T, as a function of temperature T, for various sys- 
tem sizes. For low temperatures (T < lO""*), IM2I ap- 
proaches the value |Af2| ^ 0.5 and never decays, even in 
extremely long simulations (times as long as 10® in proper 
unit^. On the contrary, in the high temperature regime 
(T > 10"^) the bi-cluster never forms, \M2\ ^ for all 
times and the PDF remains flat, as in Monte-Carlo simu- 
lations. In the intermediate regime, 10^"* < T < 10^^, the 

^ The appropriate linear time scale of the system is 2-k = 
0(1) 
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Fig. 1. (a) (resp. (b)) shows the instantaneous value of \M2\ 
as a function of time at early (resp. later) times, (c) (resp. (d)) 
shows the evolution in grey scales of the distribution 'P{6, t) at 
early (resp. later) times. The 9 axis is horizontal, the time axis 
vertical pointing downwards. 



bi-cluster forms but is less well-defined, corresponding to 
progressively smaller values of \M2\. As the temperature is 
increased there is a smooth transition to the homogeneous 
state observed in the canonical ensemble. In the next Sec- 
tion we will further discuss the internal structure of the 
bi-cluster. 

The choice of initial conditions discussed above is dic- 
tated by the need to progressively increase the energy 
starting from a ground state. One may wonder what hap- 
pens for more general initial conditions. We therefore looked 
at statistics of the final configuration, starting from many 
different realizations of a uniformly random distribution 
of angles and zero momenta. For N = 100, at temper- 
atures around T ~ 10~^ — 10"'*, the bimodal state was 
always reached over 500 initial states. For N = 1000 and 
T ~ 10~^, the totality of 250 random initial states went 
to the bimodal state. It means that the bi-cluster is fully 
attractive for this class of initial states, which is quite 
generic. 



As the formation of the bi-cluster is a unique property 
of the deterministic microcanonical system, it is an ex- 
ample of violation of the equivalence hypothesis between 
microcanonical and canonical behaviour. One might there- 
fore expect this inequivalence to show up in other mea- 
surable quantities. This is indeed the case, for example we 
find that the standard relationship between energy and 
temperature is modified. In the canonical ensemble, as 
there are only two harmonic modes, equipartition of en- 
ergy at low temperature leads to the energy temperature 
relation E = NT/2+T {N quadratic modes for the kinetic 
energy and 2 for the potential energy), which implies that 
T = 2U + O(iV-i), see (@). The eventual Order by Disor- 
der averaging over parts of the phase space with only one 
(j^adratic modes (or any finite number of them) produces 
(Sily 1/N corrections to this relation. In the microcanon- 
ifcal simulation, in the presence of the bi-cluster we find 
the anomalous relation T ~ 1.3U as shown in Fig. §(b). 
Once the temperature exceeds T > 10~^, the bi-cluster 
does not form and the energy-temperature relation crosses 
over to the canonical expression. In the presence of the bi- 
cluster, the kinetic energy is therefore much smaller than 
one should expect and the mean potential energy far in 
excess of that predicted by equipartition for the single 
non-zero mode associated with the bicluster. In fact, from 
this result, we find that the configurational contribution 
to the heat capacity d k<V > /T ^ (0.357V)/1.3 is an 
extensive quantity, in complete contradiction to the pre- 
dictions of mode counting in the canonical ensemble. 

The system manages to put a macroscopic amount of 
energy in the single non-zero mode corresponding to the 
vibrating motion of the two groups of rotators and gives 
the system a low-dimensional aspect with the bi-cluster 
taking on many characteristics of a two-particle system. 
Fluctuations in the magnetization \Mi \ , whose mean value 
is identically zero in the ground state illustrate this point. 
In zero field the potential energy is = A^|Mip/2, hence 
the extensive nature of the configurational heat capacity 
implies that mean value of the magnetization < \Mi\ > 
should be of order unity and not of order as one 

might expect in an uncorrelated paramagnetic system. 
We have indeed observed ordered features in the time 
dependence of Mi which are consistent with the non- 
fluctuational value of this quantity. 

5 Structural and Dynamical Details of the 
Bi-Cluster 

The bi-cluster never becomes perfectly formed, with \M2 \ = 
1, even in the limit [/ — > 0. Rather, the angular distribu- 
tion of rotators retains a width, with a certain population 
of rotators homogeneously distributed. In Fig. p| we show 
the angular PDF V{6) for a system of iV = 10 particles 
and three different temperatures: one very low (Fig. ^(a)) 
and two in the transition region (Fig. ||b) and (c)). The 
distribution is peaked for two angles, separated by tt, but 
there remains non-zero probability of finding angles be- 
tween the two peaks. The degradation of the bicluster 
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is progressive as the temperature is increased. We have 
made the experimental observation that the moments of 
the PDF assume the following values as the energy is de- 
creased 

< \Mk\ > ~0 for odd k 

< |Mfe| > = l/\k\ for even k with < |Mo| 1 ,(10) 

as shown in the insets of Fig. || for the first k = 20 modes. 
This observation offers no contradiction to our finding that 

< \Mi\ > is independent of system size. As the bi-cluster 
is stable at low temperatures only, the numerical value of 

< I Ml I > remains small, even though it is an intensive 
quantity. By summing the Fourier series for V{9), 
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Fig. 2. (a) \M2 \ as a function of the temperature T for N — 200 
(triangles), = 500 (circles), iV = 10^ (diamonds), iV = lO" 
(squares), (b) T vs. U for the same values of A''. The solid line 
is the modified microcanonical relation T — 1.3 U, whereas 
the dashed one corresponds to the canonical relation T = 2 U 
and is valid in the high energy regime, where the cluster has 
disappeared {{Ah] — 0). 



in Fig. H(a), where the state at time t 
iV = and T = f.43 • 10"^ is displ; 



This analytical formula is superimposed on to the numer- 
ical data in Fig. ||(a) with no free parameter apart from a 
shift of 9, due to the motion of the bi-cluster (see below). 
The agreement is impressive, and although we have no 
theory for this result, we may well say we have a solution. 
We expect the analytical formula (12) to become exact in 
the [/ -> hmit. 

The internal structure of the bi-cluster is further shown 

2.5 • 10^ with 
splayed in the {9,p) 
single-particle phase-space (so-called Boltzmann /i-space) . 
A large fraction of particles get stuck in the bi-cluster 
and have small values of \9\, while others develop much 
larger 9 values. This illustrates why we speak of a "gas- 
cluster" coexistence: particles appear to belong to two dis- 
tinct groups; those in the bi-cluster, which perform oscil- 
lations around the two centers separated by angle tt, and 
those in the gas, which have ballistic dynamics and travel 
over large distances. Once projected onto the p-axis the 
distribution is symmetric (giving, for instance, zero aver- 
age momentum), but as we see it in the {9,p) space it is 
skewed, with, on average, bigger momenta for those par- 
ticles which have traveled the furthest. 

This picture is further clarified if we fold the 0-axis 
onto [0, 27r] (Fig. |^b). The oscillating motion of the parti- 
cles of the cluster is now revealed by the spiral arms visible 
in the bi-cluster centers (see the inset). The particles in 
the cluster appear now to be distributed in a sinusoidal 
band of lesser density. Indeed, if we look inside the region 
around 9 = (Fig. m), we see that the particles are quite 
well distributed along the phase line 



p = A sin( 



(13) 



where </) is a sliding phase related to the bi-cluster motion 
and A an energy related amplitude. This gives us a hint 
as to how to discover the dynamical mechanism at work 
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Fig. 3. ViO) for N = lO" and three different values of T: 



1.46-10-^ (a), 2.6-10" 



2.8 • 10"-* (c). The full line in (a) 



is the analytical formula (112), the diamonds are the numerical 
histograms (in (b) we join them for convenience with a dashed 
line). In the insets we show the corresponding momenta of 
the PDF (diamonds), the full line is 1/n. In (a) the moments 
respect the 1/n law for even n, while in (c) only the mode 
n = 2 is non-zero. Odd modes are much smaller. 



Fig. 4. Single particle phase-space for iV = 10'', U = 1.13-10"^ 
and t = 2.5 • 10^ In (b) the angle 6 is folded onto [0, 2tv[ and 
the inset emphasizes the structure of the spirals arm around 
one of the bi-cluster center. 



in the process of particle evaporation from the bi-cluster. 
In fact, differentiating ([l^ ) with respect to time and using 
^ p, one gets for u — 29 the equation of motion 



A sinw 



(14) 



which describes a pendulum with gravity pointing up- 
wards, u — being a saddle-point. This system is inte- 
grable and cannot produce the phase-distribution in Fig. ^(a). 
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We must, therefore introduce some sort of perturbation 
capable of producing the sinusoidal layer observed in Fig. ^(b) . 
A well known mechanism for producing a stochastic layer 
in pendulum motion is the introduction of a finite time 
step, as done in the Chirikov standard map. We have 
therefore decided to iterate the map 



0.015 



u' = u + At pu 
p[, = Pu + A^At u' 



(15) 



with At — 1, taking an ensemble of initial points homoge- 
neously distributed in a small square around the unstable 
point M = = 0. A snapshot of these points at itera- 
tion 3000 is shown in Fig. ||(b); the particles have been 
stretched along the separatrix layer and distributed along 
it in a similar way to that shown in Fig. ^(a). Two features 
are responsible for the stretching: the unstable manifold 
distributes the particles away from the saddles, while the 
stable manifold attracts them to the saddles. The inho- 
mogeneous distribution of the points is a result of these 
two mechanisms. One feature of the evaporation in the full 
model, which is not at all reproduced by the map (|l5|), is 
the fact that time evolution seems to select only one of 
the two lobes of the separatrix; in other words, the sym- 
metry p —^ —p seems to be broken in the full model, while 
stretching in both the lobes is present in the map (15). 

This map is of course unable to explain why the bi- 
cluster forms and why it is stable. That is, if the evap- 
oration mechanism were the only process present then 
the bi-cluster would be progressively depleted. This is in 
contrast with numerical simulations for low temperature 
values, which show the temporal stability of this collec- 
tive state. Hence, a condensation mechanism on the bi- 
cluster should also exist, in order to establish a steady- 
state. Indeed, Fig. ^(b) shows that in the single-particle 
phase space two elliptic points are present at the centers 
of the bi-cluster, therefore in the bordering chaotic region 
of these two points a "trapping" effect could be present. 

If formula (|l^) is approximately verified, then the knowl- 
edge of the PDF at low temperatures (|l2|) would allow us 
to compute the relationship between the temperature and 
the constant A. Performing the integral, it turns out that 



2 A^ 
T=<P >=- 



(16) 



This relation is very well verified numerically for low tem- 
perature values. Moreover, particle momenta appear to 
follow a distribution which, although non-Gaussian, has 
a variance given by ([l^). The distribution shows a sharp 
peak in the center, due to the core of the bi-cluster, which 
is far from Gaussian. However, the wings of the distri- 
bution, which are due to the "gas" do follow a Gaussian 
distribution. 

One might propose that total momentum conservation 
plays a role in the observed phenomenology. An easy way 
to test this idea is by adding a small external magnetic 
field h, which removes the constraint. We observe that, 
for small h, the bi-cluster still forms and is stable at low 
temperatures. The conservation law is not therefore rele- 
vant for cluster formation. On increasing h, the distribu- 
tion 1^(0) is modified, as shown in Fig. 6. The bi-cluster 
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Fig. 5. (a) Zoom of the central region of Fig. ^a, showing the 
phenomenological law p = A sm{9 + (j)) ( p^ with A — 0.01 and 
4> = 0.3. (b) Phase-space of lO** initial particles in the square 
(u, Pu) = [-TT * 0.0005, TT * 0.0005]^ of the toy model (|l|) (with 
At = l) at iteration 3000. 



lies along the field axis and becomes asymmetric, with the 
number of rotators lying parallel and antiparallel to the 
field direction becoming unequal. The antiparallel cluster 
is continuously depleted until a single cluster eventually 
forms. This route towards a single cluster is somewhat 
counter intuitive if one thinks of a Neel ordered, unfrus- 
trated antiferromagnet. Such a system would minimize its 
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Fig. 6. P(&) in the presence of an external field for a chain 
of N = 200 particles. The diamonds correspond to h = 0.5 
whereas the triangles to /i = 0.1 and the squares to ft = 0.01. 



energy by aligning the bi-cluster perpendicular to the field 
and relaxing the two halves continuously out of the plane 
in a symmetric way. The bi-cluster would be destroyed by 
the two clusters, of equal size, rotating continuously onto 
the field direction. 



6 Moment dynamics 

As we have already remarked, there is no implicit con- 
tradiction between the proposed formula for the moments 
of the PDF ( |l0|) and the presence of a modified energy 
temperature relation. The latter means that < \Mf\ > is 
an intensive quantity rather than being 0{1/N) as in the 
in the canonical ensemble. Hence, although all the other 
odd modes of the PDF vanish in the N ^ oo limit, this is 
not true for < \Mi\ >. Rather, it remains finite in the low 
temperature regime, increasing linearly U. Solution ( p^ 
is therefore exact, only when U ^ 0. However, because 
of the smallness of U in this regime, < \Mi\ > is a small 
quantity with respect to < IM2I > and equation ( |l0| ) is 
perfectly valid. 

We expect that Mi will display interesting dynami- 
cal behaviour, given its rather unexpected extensive na- 
ture, when the bi-cluster is formed. The (Re(Mi),Im(Mi)) 
phase-plane is shown in Fig. ^(a) for N — 10* and U — 
1.13 • 10~^(T = 1.43 • 10"^). Successive points are joined 
by lines to show the relevant properties of the motion. The 
phase-point has a fast oscillatory motion through zero and 
a much slower rotatory motion centered on zero. The fast 
motion is due to the vibrations of the bi-cluster around 
the equilibrium positions (the two components of Mi cross 
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zero in phase). The slow rotational motion is due to the 
rigid rotation of the bi-cluster. This latter motion is fur- 
ther revealed by the dynamics of Af2 in Fig.|^b). Again we 
join successive phase-points with a line, showing that the 
phase point in the M2 plane is rotating around the center, 
maintaining a fixed radius \M2\ ~ 0.5. The time depen- 
dence of the phase ip2 of M2 is such that, over the long 
time span, the average is zero, but preliminary measure- 
ments of the variance show that the motion is ballistic, 
rather than diffusive, ^ t^. 



Dauxois, Holdsworth and Ruffo: Violation of ensemble equivalence in the antiferromagnetic mean-field XY model 



9 



This picture is not only qualitative, but also quantita- 
tive. Indeed, if we average the fast oscillatory motion of 
I Ml I in Fig. 0(a), we get an estimate of < >. Taking 

< |Mip > /U « 0.7, we get the correction to the energy 
temperature relation T « (2 — 0.7)J7 = 1.3C/ as shown in 
Fig.|(b). 

At the transition temperature where IM2I is decreas- 
ing, both the motion of Mi and that of M2 become more 
erratic, revealing the beginning of the region where the 
bi-cluster is progressively depleting in time. 

The other odd moments of the PDF always show an er- 
ratic motion around zero, with the variance of the cloud of 
points decreasing as N is increased; i.e. higher odd modes 
are not intensive. The higher even moments of the PDF 
show a pattern similar to that of Fig. 0(b), with a pro- 
gressively reduced radius; the motion of the higher phases 
has also ballistic features. 



7 Conclusions and perspectives 

The antiferromagnetic mean-field classical rotator system 
is an ideal laboratory to study the relation between mi- 
crocanonical and canonical ensembles. Although it has a 
trivial canonical solution, the randomly uniform state at 
all temperatures, the high degeneracy of its ground state 
induces nontrivial dynamical effects, which are revealed 
in the microcanonical ensemble. Instead of maintaining 
a random distribution of the rotators, the Hamiltonian 
dynamics selects a bimodal state, where the rotators are 
oriented along angles at distance tt with some spread, a 
bi-cluster. We have introduced an order parameter which 
reveals the formation of this state in the low energy phase 
and we have empirically obtained an analytical formula 
for the probability distribution function in angle, which 
perfectly fits the numerical data. 

In addition to this first remarkable difference between 
the two ensembles, we have also shown that the energy 
temperature relation is modified and, although still linear, 
it has a different slope in the two ensembles. The origin of 
this behavior lies in the extensive amount of energy which 
Hamiltonian dynamics puts into the oscillatory vibrating 
motion of the bi-cluster. 

Mean-field models are characterized by self-consistency. 
Indeed, a feature of our model is that rotators generate 
themselves the mean-field in which they move. Therefore, 
in the large N limit, the dynamics of our model should 
depend only on the interaction of the single rotator with 
the mean-field. As already claimed for other models (e.g . 
beam-plasma instability and vorticity defect model |16| ) 
self-consistency effectively reduces the number of degrees 
of freedom. This is why many properties of the dynamics 
of such a complex iV-body system as ours resemble those 
of a "simple" perturbed pendulum. This is also the ori- 
gin of the ordering of the rotators in a bi-cluster. There 
should be entropic reasons why the rotators prefer this 
state rather than choosing the disordered state, which is 
instead selected in the canonical ensemble. 

Many questions remain to be explored, but the most 
compelling one concerns the careful description and expla- 



nation of the dynamics. The perturbed pendulum analogy 
must be further investigated and a thorough study of the 
time evolution of the moments of the probability distribu- 
tion function in angle should allow a better understanding 
of the low-dimensional properties of the dynamics. 
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